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A considerable effort has recently gone into the study of Gaussianity of cosmic microwave 
background (CMB) data. Among such attempts, there is one with two non-Gaussianity 
indicators, proposed by the authors, and used in a search for significant deviation from 
Gaussianity in the WMAP internal linear combination (ILC) and in the single frequency 
WMAP maps with the KQTo mask. Here we extend and complement these results by 
performing a similar analysis for the de Oliveira-Costa and Tegmark (OT) WMAP three- 
year KQ75 masked map, in which the foreground is reduced through a different statistical 
cleaning procedure. 
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1. Introduction 

A key prediction of the so-called slow-roll inflationary models is that they can- 
not generate detectable non-Gaussianity of the cosmic microwave background 
(CMB) temperature fluctuations within the accuracy of the Wilkinson Microwave 
Anisotropy Probe (WMAP) P There are, however, a number of inflationary models 
that predict non-Gaussianity at a level detectable by the WMAP experiment (see 
the reviews Refs. [2]). Thus, non-Gaussianity in the CMB data may be potentially 
important to discriminate classes of inflationary models. 

The WMAP team have found that the CMB data are consistent with Gaussian- 
ity ^However, a number of recent analyses of CMB data performed with different 
statistical estimators have provided indications of deviations from Gaussianity in 
the CMB temperature fluctuations (for a list of references see, for example, Ref.4 
and references therein). On the other hand, one does not expect that a single sta- 
tistical estimator can be sensitive to all possible forms of non-Gaussianity that may 
be present in WMAP data. Thus, it is important to employ alternative indicators 
in order to shed light on the possible causes for the reported non-Gaussianity of 
CMB data. 
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Recently we proposed two new non-Gaussianity indicators, based on skewness 
and kurtosis of large-angle patches of CMB maps, which provide a directional mea- 
sure of deviation from Gaussianity on large angular scales.EI A distinctive feature of 
these indicators is that they make it possible to construct sky maps of Gaussianity 
from the CMB temperature data. In Ref. [4] we also built and studied in details such 
maps of Gaussianity generated from the data of WMAP three and five-year internal 
linear combination (ILC) maps with KQ75 mask, and found no statistically signif- 
icant deviation from GaussianityQ Here we extend and complement these studies 
by using our indicators to carry out a new analysis of Gaussianity with the CMB 
data of the so-called OT map^ three- year KQ75 masked map, in which the Galactic 
foregrounds are cleaned with the TOH procedure!^ 

2. Indicators and Maps 

The steps in the construction of our non-Gaussianity indicators and the associated 
maps, are the following^ 

i. Take a discrete set of points {j = 1, . . . , N c } homogeneously distributed on the 

CMB celestial sphere S 2 as the centers of spherical caps of a given aperture 7; 
ii. Calculate for each cap j the skewness and kurtosis given, respectively, by 



where N p is the number of pixels in the j cap, is the temperature at the i 
pixel, T is the CMB mean temperature, and a is the standard deviation. Clearly, 
the numbers Sj and Kj obtained in this way for each cap can be viewed as a 
measure of non-Gaussianity in the direction of the center of the cap (9j, cf>j). 
iii. Patching together the Sj and Kj values for each cap, one obtains our indicators, 
i.e., discrete functions S = S(9,(f>) and K = K(9,<f>) defined over the celestial 
sphere, which can be used to measure the deviation from Gaussianity as a func- 
tion of (9, 4>). 

Clearly, the discrete functions S = S(9, <f>) and K = K(9, <j>) can be expanded 
into their spherical harmonics in order to have their power spectra S( and Kg. 
Thus, for example, for the skewness one has S{9, <f>) = X^fco ^2m=-£ ^ m Ytm(&i 4>) 
and Si = (2£ + l)" 1 J2m \^im\ 2 - Similar expressions obviously hold for the kurtosis 



a We have also used our indicators to test for Gaussianity the single frequency foreground unre- 
moved WMAP three and five-year maps, but here for comparison we focus on the foreground- 
reduced ILC and OT maps. 

b For a detailed discussion of the indicator briefly presented here we refer the readers to Ref. [4] 




(1) 



K = K{9,4>)- 
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3. Main Results and Conclusions 

In this section we shall report the results of our calculation of S = S(8, <j>) and 
K = K(6, (j)) indicators for the WMAP three-year data of the OT with KQ75 
mask, and compare with the results of Ref. [4| for the ILC WMAP three- year map 
with the same mask. 




Fig. 1. Monopole removed Skewness (left) and monopole removed Kurtosis (right) maps gener- 
ated from the WMAP three-year OT map with mask KQ75. Colored high resolution versions of 
these figures are available in the arXiv version. 

In order to minimize the statistical noise, in the calculations of skewness and 
kurtosis indicator maps (S— map and K— map) we have scanned the celestial sphere 
with spherical caps of aperture 7 = 90°, centered at N c = 12 288 points on the 
sphere homogeneously generated by using HEALPix package.^ 

Figure [T] shows the Mollweide projection of the monopole removed S (left) and 
K (right) maps generated from the OT three-year map with a KQ75 mask. These 
maps show spots with higher and lower values of indicators S{9,4>) and K(9,<f), 
thus suggesting at first sight non-Gaussianity of the CMB as given by OT tem- 
perature map. Fig. [1] is also suggestive of large-scale dominant components (low i 
components) in the maps of both indicators. 



S, : ILC-3yr & OT-3yr 
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K, : ILC-3yr & OT-5yr 
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Fig. 2. Low I differential power spectra of skewness \Si — Sf\ (left) and kurtosis (right) \Ki — K(\ 
calculated from the monopole removed maps of Fig. [T] which were generated from the WMAP 
three-year OT map with mask KQ75. The 68% and 95% confidence levels (obtained from MC 
simulated radom maps) are indicated, respectively, by the dashed and dash-dotted lines. 



In order to obtain quantitative information about the observed inhomogeneous 
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distribution for these non-Gaussianity indicators, we have calculated the power 
spectrum of both the S— map and K— map for the three-year OT KQ75 cut-sky 
map. These power spectra allow to estimate the statistical significance of the S 
and K multipole values, by comparing their power spectra with the averaged power 
spectrum of the S and K maps obtained by averaging over a set of 1 000 Monte- 
Carlo-generated statistically Gaussian CMB maps. In the left panel of Fig. [2] we 
depict differential power spectrum \Se — Se\ (over line denotes the mean) are obtained 
from the Monte-Carlo simulations with each map being a stochastic realization 
of WMAP best-fitting angular power spectrum of the ACDM model, obtained by 
randomizing the multipole temperature components ai m within the cosmic variance 
limits. The right panel shows the corresponding plot \Kg — Ki\ for the K— map. For a 
comparison the panels of Fig. [5] also show the differential power spectrum obtained 
from the WMAP three-year ILC KQ75 cut-sky map. 

Figure [5] shows the deviation of St and Kg values (£ = 1,- • - ,10) from the mean 
values of Si and Kg are within (less than) 95% the mean values (for each £) ob- 
tained from 1 000 Monte-Carlo-generated statistically Gaussian maps. This clearly 
indicates how likely is the occurrence of the multipole values obtained from the 
OT data in the set of MC multipole values (from MC maps), giving therefore a 
clear indication of consistency with Gaussianity for the OT KQ75 masked map, in 
agreement with previous analyses performed with different statistical tools (see, e.g. 
Refs.HSl). 
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